On the influence of noise on chaos in nearly Hamiltonian systems 
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The simultaneous influence of small damping and white noise on Hamiltonian systems with chaotic 
motion is studied on the model of periodically kicked rotor. In the region of parameters where damp- 
ing alone turns the motion into regular, the level of noise that can restore the chaos is studied. This 
restoration is created by two mechanisms: by fluctuation induced transfer of the phase trajectory 
to domains of local instability, that can be described by the averaging of the local instability index, 
and by destabilization of motion within the islands of stability by fluctuation induced parametric 
modulation of the stability matrix, that can be described by the methods developed in the theory 
of Anderson localization in one-dimensional systems. 

PACS numbers: 05.45-a, 05.40-a 



I. INTRODUCTION 

From the point of view of chaotic dynamics, the Hamil- 
tonian systems are marked out by the omnipresence of 
chaos: for nearly any Hamiltonian system with not less 
than one and a half degrees of freedom (with the exemp- 
tion of completely integrable models that are non-robust 
and therefore exceptionally rare) the chaotic motion is 
possible for some initial conditions. On the contrary, 
for the dissipative systems of the same complexity of 
the structure chaotic motion on strange attractors cither 
could be attained only in limited domains of the param- 
eter space or is inaccessible at all 0, • 

Inclusion of the dissipative terms, even arbitrarily 
small, in the canonical equations of motion of the Hamil- 
tonian system can change the character of the motion 
drastically. In particular, such addition can banish the 
chaos: for example, for the autonomous Hamiltonian sys- 
tems with added (viscous) damping the only possible at- 
tractors are stable fixed points. It must be noted that 
this abrupt change may be basically formal, resulting 
from the presence of the transition to the infinite time 
limit in the rigorous definitions of important character- 
istics of chaotic motion, like the Lyapunov exponent and 
correlators of dynamic variables. In many experimen- 
tally relevant models the ratio of the dissipation 7 to the 
typical frequency of motion lu may take very small val- 
ues. Thus, for radiation damping of vibrations of poly- 
atomic molecules one has j/oj ~ 10~ 10 ; the same order 
of magnitude of j/lo turns out for the tidal friction of the 
celestial bodies of the Solar system. In these situations 
the duration of the "transient chaos" phase T ~ 7 _1 is so 
long that accurate determination of characteristics of the 
chaotic motion can be carried out without the account of 
dissipation. 

Physically the introduction of dissipation in the equa- 
tions of motion is a form of description of the interac- 
tion of an isolated (in the zeroth approximation) system 



'Electronic address: pve@shg.phys.msu.su 



with its environment - a "heat bath" with practically 
infinite number of degrees of freedom, continuous spec- 
trum of eigenfrequencies and internal dynamics that is 
independent of the state of the selected system. This 
heat bath may be considered also as a source of noise 
- that is, acting on the selected system random forces, 
whose statistical characteristics are determined by the 
properties of the heat bath alone. The problem of simul- 
taneous influence of small dissipation and weak noise on 
the features of chaotic motion in the originally Hamilto- 
nian non-autonomous system is the main concern of this 
paper. 

The studies of the influence of noise on chaotic motion 
were pioneered by Lieberman and Lichtenberg Q just 
by analysis of the effect of fluctuations on the Hamil- 
tonian non-autonomous system. However, the modern 
paradigm of the domain was formed later by Crutchficld 
and Huberman 0, who switched the attention to the 
exploration of strongly dissipative systems (sec late re- 
view |6()- The influence of noise on the Hamiltonian sys- 
tems has been discussed recently in the context of the 
problem of decay of metastable chaotic states 0, , but 
in general the field doesn't seem to be fully investigated. 

On the contrary, the influence of small dissipation on 
the Hamiltonian chaos is well understood: Afraimovich, 
Rabinovich, and Ugodnikov £ 9] have shown, that with 
switching on a small dissipation phase trajectories of sta- 
ble periodic motions of the Hamiltonian non-autonomous 
systems become attractors with regular motion, and 
chaos disappears. With the further increase of 7 these 
attractors may lose their stability; annihilation of the 
last one turns the system back into chaotic motion on a 
strange attractor, that resembles the chaotic motion of 
the original Hamiltonian system. This pattern needs two 
specifications. First, the strange attractor may emerge 
before vanishing of the last of regular ones - the system 
could be multistable. This case, mentioned in Q as "log- 
ically possible" , will be met in our model. Secondly, if 
the Hamiltonian system has no islands of stability that 
correspond to periodic motion, then the transition from 
Hamiltonian to dissipative chaos can occur immediately. 
This case, apparently, will be present in our model too. 
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The aforesaid permits to specify the main problem of 
our paper: what intensity of noise is necessary to restore 
the chaos, repressed by dissipation? 

The rest of the text is organized as follows. In Sec. 
II the basic model is introduced. Sections III and IV 
treat two mechanisms of restoration of chaos by noise: 
fluctuation transfer to domains of local instability and 
parametric destabilization of motion within stability is- 
lands. Sec. V treats the influence of strong noise on the 
Lyapunov exponent and correlation functions of chaotic 
motion. Sec. VI contains the summary of results and 
their discussion. 



II. THE BASIC MODEL 

We start from the well-known periodically kicked rotor 
the non- autonomous model with the Hamiltonian 



t2 °° 

H(I,0,t) = — +Kcos9 S(t-n), 



(1) 



where I and 9 are the dynamic variables (canonically 
conjugated momentum - action and coordinate - angle), 
K is the control parameter, and 5{z) is the Dirac delta- 
function. The stroboscopic mapping that links values 
of dynamic variables at the moments of time n — and 
n + 1 — 0, preceding two consequent kicks, 



I' = I + Ksm( 



9' =9 + 1', 



(2) 



is known as the standard, or Chirikov - Taylor, mapping 
and thoroughly studied ■ 

The generalization of the model that includes dis- 
sipation and noise will be described by the equation of 
motion for the angular variable 



9 + j9 - K sin 9 ^ 6 (* 



at), 



(3) 



n— — oo 



where 7 is the damping constant. The £(t) in the RHS 
is the Langevin random force, that is a stationary, dis- 
tributed by a Gaussian law, 5-correlated random process 
(white noise) with zero mean and the correlator 



(£ (t)£ (t + r)>= 2705 (r), 



(4) 



where = ksT is the noise temperature in energy units 
(in the system of scales of the model). The model given 
by Eqs. 10 and (@J has three parameters: K, 7 and 
0. We shall restrict ourselves by the domain of small 
damping, 7 <C 1, where the system is nearly Hamiltonian. 

In the absence of noise, at = 0, the stroboscopic 
mapping for this model is given by equations 

I' = a (I + K sin 9) , 0' = 9 + b (I + K sin 9) , (5) 

where 



a = e 



b = 1 -\l-e-'<). 



(6) 



The two-parameter mapping given by Eqs. JSJ is a spe- 
cial case of the four-parameter Zaslavsky mappi ng^ that 
has been introduced in 0] and studied in [TJ El El • 
The main efforts of these studies were applied to the case 
7 > 1. Here we describe in brief the properties of our 
model for the case of small damping, 7 <C 1. 

At small and moderate values of the control parameter 
K the most important attractor of the mapping Eqs. (J5J) 
is the fixed point I = 0,9 = %, that is stable in the range 



K<K X = 



2(1 



1 H Y 

12 ' 



(7) 



For K > Ki the leading attractor is the symmetric 
cycle C| of two points that are related by equations 
V = -I, 9' = 2tt - 9. It is stable while 



K < Ko — 



7T (1 + a) 



2,t( 1 + ^7 



(8) 



For K > Ki the leading attractors are two asymmetric 
cycles of length two C% 1 and Cf 2 . The phase coordinates 
of their points are related by equations I' = —I, 9' = 
it + 9. They are stable in the domain 



2a/tt 2 + 1 (1 + 0.1527) 



(9) 



In general case the system defined by Eqs. (JSj is mul- 
tistable. From the first of these equations it follows that 
the strip 

I'l < 1+ = — 7 (10) 

exp 7 — 1 

is the trap (the absorbing set) of the system: any phase 
trajectory that comes within this strip never leaves it. To 
determine the comparative roles of basins of attraction 
of strange and regular attractors of the model the frac- 
tion / of chaotic trajectories among the set with random 
initial conditions, uniformly distributed within the trap, 
was calculated. The results are presented in Fig. 1. It 
can be seen that the strange attractor is born within the 
domain of stability of the cycle Cf , and after the loss of 
stability of cycles Cf 1 and Cf 2 , in the range K > ^3(7), 
it remains the only apparent attractor of the system. 

For > the stroboscopic mapping for the system 
Eq. © has the form 



I' = a (I + K sin 9) + v, 
9' = 9 + b{I + Ksm9) + tp, 



(11) 



where v and if are the random increments (fluctuations) 
of action and angle for the unit interval of time. Fluctu- 
ations v and ip at the moment of time t after the begin- 
ning of the motion with definite initial conditions have 
the Gaussian distributions with the dispersions 



0(1 




-2 7t) ; 



{2~ft - 3 + 4e- 



■ft 



-2 7 n 



(12) 
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K < 5.4, when the motion of the system in the absence 
of noise is regular, the Lyapunov exponent increases with 
6 and at some value of 6o passes through zero. 

We turn to the theoretical description of the onset of 
chaos. For conservative (area-preserving) mappings with 
strong chaos rather accurate estimate for the Lyapunov 
exponent could be obtained by averaging of the local in- 
stability index - the logarithm of the maximal in abso- 
lute value eigenvalue of the stability matrix - over the 
domain of chaotic motion 0, Since our model is 

nearly Hamiltonian, we may try to use this approach. 

The stability matrix for the mapping given by Eqs. 
CU is 



FIG. 1: The dependence of the fraction / of the trap area, 
covered by the basin of attraction of the strange attractor 
of the system Eqs. ©, on the control parameter K. The 
damping value 7 = 0.05. For each value of K 100 uniformly 
distributed initial conditions was taken. 



A = 



a aK cos 
b 1 + bKcost 



(16) 



The local instability index depends only on the angle 0; 
for K < 4 



respectively [15]. In our case jt = 7 <C 1 and the expres- 
sions Eqs. (|T2)l could be replaced by their asymptotics 



D v = 2 7 9, 



(13) 



Fluctuations v and tp are positively correlated; that 
has clear physical meaning: positive increment in veloc- 
ity (that is numerically equal to the action) at the unit 
interval of time leads, most probably, to the positive in- 
crement of coordinate. The correlator M = (vip) by the 
moment t after the beginning of motion can be calculated 
by the method described in [15]: 



a (0) = In 




(17) 



where S = a + 1 + bK cos is the trace of A and D = a 
is its determinant. For small 7 almost everywhere in the 
interval ir/2 < < 3tt/2 the index is negative and con- 
stant, a{&) = —7/2, and the motion is locally stable. 
Most probably the system stays in this domain, but un- 
der the influence of noise it can sporadically enter the 
domains of local instability. For large enough values of 
their contribution can compensate the weak squeezing 
of phase trajectories in the central part of the attractor. 



M 



6 

1 



(1 



- 7 t\ 



(14) 



For jt = 7 <C 1 the asymptotic value of this correlator is 
M = O7. The joint distribution of fluctuations of action 
and angle has the form 



<(v,<p) = 



V3 
27T79 



exp 



(15) 



In the presence of noise the phase trajectory can reach 
any point of the phase space of the system. However, for 
finite 7 the system with overwhelming probability will 
stay in the strip with limited action values that is much 
narrower than the trap given by Eq. IjlUI) . For the de- 
scription of this domain of concentration of the probabil- 
ity density the term "attractor" will be used. 



III. THE THRESHOLD OF CHAOS: TRANSFER 
TO THE DOMAIN OF LOCAL INSTABILITY 

The condition of existence of chaos is, by definition, the 
positive value of the Lyapunov exponent a. Numerical 
calculation show, that in our model at moderate values of 



A. Small K 

For small values K <C 1 in the absence of damping (7 = 
0) the evolution of the periodically kicked rotor nearly 
everywhere in the phase space can be described by the 
time averaged (and thus time-independent) Hamiltonian 
of the system, that is given by Eq. : 



H{I,0) = H(I, 



' 2 



K cos 0. 



(18) 



For an autonomous Hamiltonian system the inclusion of 
viscous damping and connection to the Langevin (white 
noise) heat bath lead to the canonical distribution of 
probability in the phase space 



W (I, 0) = TV exp _ 



H(I, 

e 



(19) 



where N is the normalization constant. For the averaging 
a = (a (0)) one needs to know the angular distribution 
W{&). Its normalized form could be found from Eqs. i|18|) 
and fTCty: 



W(6) 



1 



K 



(20) 
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FIG. 2: The dependence of the temperature threshold of 
chaos 0o on damping 7 for K = 0.3. Calculation by Eqs. 
(121 1 and 1221 (line) and numerical calculation (points). 



FIG. 3: Numerically found dependence of the temperature 
threshold of chaos Go on the control parameter K for the 
damping value 7 = 0.1. 



where Iq(z) is the zeroth order modified Bessel function 
of the first kind. Since nearly all probability density is 
concentrated in the stability interval, where cos 9 < 0, 
the contribution of this domain to the averaged value is 



7 
2' 



(21) 



To calculate the contributions of zones of local insta- 
bility we neglect the damping; then we have a (9) 

cos 9. Positive contribution of two instability strips 
could be estimated by the integral 



1 



2ir 

- J cxp ^— — cos 9^j V K cos 9 d6. 

3tt/2 



tt/o (Aye) 



(22) 

The threshold values Bo found by averaging of the local 
instability index and by direct numerical calculation are 
compared in Fig. 2. 

If A/O 3> 1, the integral can be calculated analyti- 
cally: replacing the Bessel function by its asymptotics 
for large value of the argument, and approximating the 
cosine by the linear function, we obtain: 



1 n K 

2K B 6XP "0- 



(23) 



From the condition (<x) = er_ + a + = the threshold 
value of temperature is determined by the root of the 
equation 



/ 2 n K 
7 =^-e exp --. 



(24) 



This expression yields the asymptotic dependence of the 
temperature threshold of chaos for small K: it has the 
form 



0c 



A 



|ln 7 | 



(25) 



and possesses a logarithmic accuracy. 



B. Large K 

The results of the previous subsection shows that for 
small K the threshold 0o grows with the increase of the 
control parameter K, 0o oc K. On the other hand, as it 
was noted in the Sec. II (see Fig. 1), for K > K3 chaos 
in the dissipative system exists (apparently for any initial 
conditions) even without any noise. For the reasons of 
continuity we may expect that there may exist a range 
of values of A where the dependence ©o(A) at the con- 
stant damping 7 is decreasing. The numerical calculation 
confirms this suggestion (see Fig. 3). 

For the area-preserving mapping Eq. (2), reduced to 
the basic square (0 < I < 2tt) x (0 < 9 < 2tt), for large A 
the chaotic motion takes place in the chaotic component 
that covers the largest part of the phase space (for K > 2 
the measure of the chaotic component /x(A) > 0.78). For 
the system with damping and noise Eq. we shall re- 
tain the name "chaotic component" for the part of the 
attractor that includes the points of the chaotic compo- 
nent of the conservative system, and the complementary 
part will be referred to as " islands of stability" . Limiting 
ourselves to the case A < 4, we will take into account 
only one island of stability that surrounds the stable fixed 
point I = 0, 9 = 7t. 

Some properties of chaotic motion of the system with 
damping and noise in the chaotic component could be 
described by the following simple model. The action 
variable for one time step receives the increment A7 rs 
Asin# with the averaged square value (A/ 2 ) w K 2 /2. 
For the motion in the chaotic component the correla- 
tions of the consequent values of 9 are small 0, Il8| . 
and we can depict the evolution of the system as the 
motion of the rotor with the damping 7 under the influ- 
ence of some Langevin force, a white noise coming from 
the source with an effective temperature 0*. From Eq. 
(fH)l we have the estimate 0* w K 2 /Aj > 1. In this 
approximation the distribution of phase density in the 
chaotic component will become canonical one, with uni- 
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form distribution of angles and the Gaussian distribution 
of action, 



w c (i,e) 



2n 3 K 2 eX P I K 2 



27 j2 



(26) 



This expression is applicable for small 7 and large K . 

Let's assume that in the island of stability the prob- 
ability density also has the distribution of the canonical 
form 



W t (I, I 



N exp 



H(I, 



e 



(27) 




4 5 6 

inverse temperature p 



where H (I, 9) is an effective Hamiltonian (a function that 
is constant on the invariant curves of the standard map- 
ping (2)), and N is the normalization constant. This 
assumption is plausible in view of Eq. (|19[1 : additional 
support for it will be obtained in the next section. 

If 7 and are sufficiently small, then the phase tra- 
jectory can leave the island of stability or return to it 
only by passing through the narrow strip of the width 
S ~ \Ay0 along the border of the island of stability. The 
probability P of finding a phase point in the chaotic com- 
ponent could be found from the balance considerations 
by equalizing W c and Wi on this border. For P< 1 we 
can neglect the non-uniformity of the distribution Eq. 
(|26|l in action and obtain the estimate 



P = 



7 



A 

exp,-- 



(28) 



where A is the value of the effective Hamiltonian H (I, 9) 
on the border of the island of stability. At present we can 
not calculate this quantity analytically, but from its value 
taken from the numerical calculations (and depending 
only on K) by Ea. (|28|) we can find the dependence of P 
on and 7. 

In the numerical experiment the basic square 
(0 < I < 2tt) x (0 < 9 < 2n) has been separated into 10 4 
cells. A chaotic trajectory of the standard mapping Eq. 
(J2J) in this square was calculated for 10 5 time steps, and 
all cells, in which the trajectory came at least once, were 
marked as the mask of the chaotic component. Then for 
the trajectories of the system with damping and noise 
Eq. the probability P has been calculated as the frac- 
tion of points of the trajectory whose projections on the 
basic square got into one of the cells of the mask. 

The results of numerical calculations for the values 
K = 3 and 7 = 0.05 are shown in Fig. 4. Fit of the linear 
dependence of InP on the inverse temperature (3 = _1 
for these points gives values A = 1.07 and N — 0.057. 

From the assumption that the motion inside the sta- 
bility island gives to the Lyapunov exponent the nega- 
tive contribution cr_ = —7/2, and that the positive con- 
tribution from the motion in the chaotic component is 
cr_l_ = a (K) P, where <j{K) is the Lyapunov exponent 



FIG. 4: Dependence of the probability P of finding a trajec- 
tory of motion of the system with damping and noise in the 
place of the chaotic component of the standard mapping on 
the inverse temperature (3 = B _1 for the values K — 3 and 
7 = 0.05. Numerical calculation (points) and linear fit to the 
points (line). 



value of the Hamiltonian system, for the temperature 
threshold of chaos we obtain the equation 



0o 



A (J?) 



In 



a {K) N^J^K 2 ^ 



(29) 



The threshold values 0o found by calculating the prob- 
ability of transfer to the chaotic component by this for- 
mula and by direct numerical calculation are compared 
in Fig. 5. From the Eq. Ij29(l it is seen that essentially 
©o is proportional to the "activation energy" A(K); the 
dependence on other parameters is only logarithmic. The 
general behaviour of the dependence on Fig. 3 could be 
explained by decrease of the size of the stability island 
with the increase of K; the dip around K — 4 reflects the 
restructurization of the regular attractor. 




0.00 



0.02 



0.04 0.06 
damping y 



0.08 



0.10 



FIG. 5: The dependence of the temperature threshold of 
chaos 0o on damping 7 for K — 3. Calculation by Eq. I29H 
(line) and numerical calculation (points). 
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IV. THE THRESHOLD OF CHAOS: 
PARAMETRIC DESTABILIZATION IN THE 
ISLANDS OF STABILITY 

Although the agreement between the theoretical curve 
and the numerical points in Fig. 5 is rather convincing, 
the increase of discrepancy at very small 7 is strange: 
just in this domain the damping must have especially 
little influence, and the picture of transfer between the 
island of stability and the chaotic component promises 
to be asymptotically exact. Furthermore, this discrepan- 
cies could not be neglected from the quantitative point of 
view, since for the sharp dependence of -P(0) small vari- 
ations of O at low temperatures produce large changes 
in the positive contribution to the Lyapunov exponent 
(T + = a (K) P. For example, for K — 3 and 7 = 0.01 
substitution of the numerically found value 0o = 0.086 
in the Eq. gives P = 2.6 -lO" 5 and a + = 1.8 -10" 5 = 
3.5 • 10~ 3 (7/2). Thus we must conclude that there exists 
another mechanism creating the instability that acts on 
the parts of phase trajectories that are localized within 
the islands of stability. 

In the domain 2 < K < 4, in which the island of 
stability surrounds the fixed stable point / = 0, 6 = ir, 
we can consider the dynamic variables x = 8 — tt and 
y = I to be small. Then by substitution cos x « 1 — x 2 /2 
the matrix of stability could be represented in the form 



.4 



a —aK + arj 
b 1-bK + brj 



(30) 



where 77 — Kx 2 /2 are small corrections. The quantities 

77 are fluctuating under the influence of noise and could 
be treated as random. The stochastic modulation of pa- 
rameters of the mapping leads to destabilization of the 
motion. The corresponding exponent of instability + 
we will calculate in the conservative approximation, since 
the contribution to the common Lyapunov exponent from 
damping er_ = —7/2 and from stochastic modulation cr+ 
for small 7 are additive. 

The transformation of variables determined by matri- 
ces Eq. (|3"U)l at 7 = can be reduced to the three-term 
recurrent relation for the angular variable: 



x n+ i - (2 - K - Tj„) x n + X n -i = 0. 



(31) 



This expression can be interpreted as an equation for am- 
plitudes x n of the stationary wave function in the quan- 
tum one-dimensional tight-binding model with unit non- 
diagonal matrix elements (transfer integrals) between ad- 
jacent sites, random site energies r\ n and the energy eigen- 
value E — 2 — K (one-dimensional chain with diagonal 
disorder). The calculation of the Lyapunov exponent for 
this system was carried out in the context of the theory 
of Anderson localization. For the case in which r\ n are 
independent random variables with zero mean, (77) = 0, 
and small dispersion, \(f) -C 1, the Lyapunov exponent 



has been calculated by Derrida and Gardner 0] . Corre- 
lations between consequent values of r\ n were taken into 
account by Tessieri and Izrailev [2£j ■ The stochastic Lya- 
punov exponent is proportional to the dispersion of fluc- 
tuating parameter 77: 



2 {AK - K 2 ) 



C(w). 



(32) 



The correlation factor C{lo) in this expression has the 
form 



C (u) = 1 + 2 



00 

b v (k) cos (2wfc), 

fe=i 



(33) 



where b v (k) are normalized correlation functions of the 
random variable 77, b n (k) = (rjirji+k) / ( 7 ? 2 )) an d 



K 



(34) 



is the average angle of rotation of a vector by the lin- 
earized standard mapping. Formula Eq. (|32|l is applica- 
ble for the values of K that are not too close to K = or 
K = 4. In our case (77) ^ 0, but we can include this value 
in the parameter K; this renormalization will change its 
value into K = K (l — ^£ 2 )/2). In what follows we shall 
retain the designation 77 for the fluctuating quantity with 
zero mean, 77 = K (x 2 — (a? 2 )) /2. To use the expression 
Eq. H32|) we have to determine statistical characteristics 
of the variable r\: its dispersion and correlation function. 
It must be noted at once that for calculation of these 
quantities the account of damping is essential. 



A. Invariant density in the island of stability 

Since for low temperatures O the phase trajectory 
nearly all the time is located in the vicinity of the sta- 
ble fixed point, for finding the invariant distribution of 
the probability density W{x, y) we may use the linearized 
mapping L, 



x 1 = by + (1 - bK) x, y' 



ay ■ 



aKx. 



(35) 



The motion of the system on a unit time interval can 
be separated into two stages: the first is the evolution 
under the mapping Eq. (|35|1 . the second is addition of 
the fluctuation increments (cf. Eq. I|llf> ). The probabil- 
ity of coming in the vicinity of the point (it, v) after the 
first stage is proportional to the value of the invariant 
density in the vicinity of its prototype L^ 1 (u,v). The 
influence of noise can be described by the convolution of 
the obtained distribution with the distribution of fluctua- 
tion increments w(v, ip) Eq. (|15fl . Thus for the invariant 
density we obtain the following integral equation: 
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W{x,y) = a ^ J J dudvW 



b T, 1 - bK 

—v,Ku-\ V 

a a 



_V3_ 

27T76 



exp 



1 

76 



3 (u — x) 2 — 3 (u — x) (v — y) + (v — y) 2 



m 



We will look for its solution in the form of a canon- 
ical distribution with the effective Hamiltonian, that is 
bilinear in action y and angle x: 

W (x, y) = exp ~ (Ax 2 + Bxy + Cy 2 ) . (37) 

After substitution of this expression in Eq. (|36|l . inte- 
gration and equivalizing the coefficients at the identical 
powers of dynamic variables, in the lowest order in 7 we 
obtain the parameters of the function Eq. (|37|l : 



A 



3K 
6-K' 



D 



3K 
6-K' 



C 



6-K 



(38) 



Now we can calculate the moments of dynamic variables, 
e.g. 

, 2 = 12-2A- = 6-K 

X 1 12K-3K 2 ' K y ' 12 -m ' V ' 

and the dispersion of the fluctuating quantity, 



W) = \ ( 



1 (Y2-2K 



12-3K 

that enters in the RHS of Eq. ^ 



e 2 



(40) 



B. Correlation function 

With the known form of the invariant density, the cor- 
relation function of the angular variable x could be cal- 
culated by the direct integration. For the linearized map- 
ping Eq. I|35|l the values of B x (n) = (xjXj_|_„) could be 
expressed through two moments, (x 2 ') and (xy). For ex- 
ample, 

B x (0) = (x 2 ) B x (1) = (1 - bK) (x 2 ) + b (xy) , 

B x (2) = (l-bK(2 + a)+ b 2 K 2 ) (x 2 ) 

+b (a - bK + 1) (xy) . (41) 

For small 7 the normalized correlation function b x (n) — 
B x (n)/B x (0) can be expressed in the form 

b x (n) = cos (Qn) exp ( — • (42) 

Here the tilde over to reminds that the renormalized value 
of K must be used in calculations. This formula is com- 
pared to the numerical calculations in Fig. 6. 

The normalized correlation function of the fluctuating 
variable 77 can be calculated in a similar way: 
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FIG. 6: Normalized correlation function b x (n) of the angular 
variable for the values of parameters K = 3 and 7 = 0.2. (a) 
For 9 = 0.05, (b) for 9 = 0.2. Calculation by Eq. JEJ (line) 
and numerical calculation (points). 



With this expression one can calculate the correlation 
factor C(uj) (see Eq. I|33|l '). For small 7 it is given by the 
expression C (lo) rj (27)" . With it from the condition 
cr_ + (j + =0 follows the estimate of the temperature 
threshold of chaos 



9o = 2^AK-K? 12 3i V 
V 12 - 2K ' 



(44) 



br/ (n) = cos (urn) exp (—771) . 



(43) 



This formula is too crude for the practical application: 
it gives only the estimate of 0rj from below. Here is 
the reason for this limitation: the expressions Eqs. I|42() 
and (|43[) for the correlation functions are valid only for 
small temperatures, < 7. For larger values the non- 
linear terms that are present in the exact mapping Eq. 
(|llfl change the frequency of oscillations of the correla- 
tion function b„(n) (see Fig. 6(b)); by this they spoil the 
resonance with the cosine factor under the summation 
sign in Eq. ll-i-it and considerably decrease the value of 
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FIG. 7: The dependence of the Lyapunov exponent a on the 
noise temperature O for K — 3 and 7 = 0.1 obtained from 
numerical calculation (points). Parabola in the left part (solid 
line) - calculation by Eq. 14511 with k = 1. Horizontal dashed 
line marks the limiting value (Too = 0.538. 

C(uj), down to the value about 2-3 on the threshold of 
chaos. The numerical calculation shows that around Oo 
the temperature dependence of the Lyapunov exponent 
is accurately described by the formula 

cr=-| + K6 2 (45) 

with a coefficient k about unity. From Eq. I|45|l the sim- 
ple approximation for the temperature threshold follows: 

6o « c^7, (46) 

where constant c is about unity. The fit of the law Eq. 
(|4*tJ)l to the points in Fig. 5 gives c » 0.88. 
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FIG. 8: The dependence of the limiting value of the Lyapunov 
exponent aoo on the control parameter K for 7 = 0.1. Calcu- 
lation by Eq. I|47^ (line) and numerical calculation (points). 



the band edge, namely: 

o-oo = -| + 0.229X 2/3 . (47) 

The dependence given by this equation is in good agree- 
ment with the numerical data up to the value of K ss 1 
(see Fig. 8). 

It may be noted that from Eq. I|47|l it follows that for a 
given value of damping 7 there exist a range of values of 
the control parameter K < 3.23 7 3 / 2 in which chaos could 
not be reached for any intensity of the noise. For large 
values K > 3 the limit Coo does not differ noticeably from 
the Lyapunov exponent cr(K) of the original Hamiltonian 
system. 



B. The angular correlations 



V. STRONG NOISE 

In this section we will look at the effects of noise with 
temperature much higher than the chaos threshold 60 ■ 

A. The Lyapunov exponent 

With the increase of the noise temperature the Lya- 
punov exponent increases monotonously and tends to 
some finite limit Coo (see Fig. 7). 

Since with the increase of noise the typical values of 
increments of the angle (</?) and of the action (v) grow, 
<p ~ v ~ ^/~fQ (see Eq. {EJl), for 9 > 7 _1 all cor- 
relations of dynamical variables vanish. Therefore the 
limiting value <Too is equal to the Lyapunov exponent of 
the infinite product of the matrices Eq. (|17f) with uncor- 
related values of 9, uniformely distributed in the interval 
< 9 < 2tt. For small K the value of can be calcu- 
lated from the result of for the localization length at 



In the theory of the standard mapping it is customary 
to study the angular correlations through the properties 
of the variable s = sin 9 0, . From the symmetry 
considerations it has zero mean: (s) = 0. 

Let's consider the correlation of two consequent values 
of this variable: 

B s (1) = (sin 9 sin 9') . (48) 

When the invariant density W(I, 9) is known, the calcula- 
tion of the correlation Eq. (|48|l is reduced to the two-fold 
integration. For small K the distribution can be taken 
as the canonical one with the averaged Hamiltonian Eq. 
(|18|l . With approximating the angular distribution by 
the Gaussian function, for small damping 7 we have the 
expression 



Bs (1) = 




It is compared to the numerical data in Fig. 9(a). 
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FIG. 9: Dependence of the correlation function B s (l) of the 
angular variable s — sin 8 on the noise temperature B. (a) For 
K = 0.3, 7 = 0.01. Calculation by Eq. I49H and numerical 
calculation (points), (b) For K — 3, 7 = 0.05. Calculation by 
Eq. 1501 (line) and numerical calculation (points). 



For large values of K > 2 and small temperatures the 
value of B s (\) could be calculated from the distribution 
Eq. 



B s (1) = (1 - K) (.x 2 ) + (xy) = 



12-8K + K' 
12K-3K 2 



6. (50) 



It is compared to the numerical data in Fig. 9(b). The 
numerical calculation shows in this case the decrease of 
B s (l) with 9 for sufficiently high temperatures. The rea- 
son for it is qualitatively clear. We have two competing 
contributions to B s (l): a negative one from the island 
of stability and a smaller positive one from the motion 
in the chaotic component. The increase of the temper- 
ature leads to the increase of the probability P of the 
motion within the chaotic component, that eventually 
suppresses the negative contribution. The way of accu- 
rate calculation of B s (l) for high temperatures at present 
is not known. 

Thence both for small and large values of the control 
parameter K the increase of the noise temperature (from 
zero) induces first the increase of the correlation of the 
consequent values of the angular variable s — sin 8 up to 
some maximal value, and then its decrease. 



VI. CONCLUSION 

Above we have studied the model of periodically kicked 
rotor with added damping and white noise. We expect 
that some of the established features and relations are 
typical and will hold at least qualitatively for many rep- 
resentative non-autonomous Hamiltonian systems with 
chaotic motion. With this in view, in this Section we 
summarize main results of this paper in a generalized 
way. Since the growth of the control parameter of the 
standard mapping K increases both the amount (given 
by the invariant measure fi of the chaotic component) and 
the intensity (given by Lyapunov exponent a) of chaos 
we will refer to K as to the strength of chaos. In what 
follows (as well as everywhere above) the Lyapunov expo- 
nent denotes the largest of the characteristic exponents 
of the stroboscopic mapping of the system that can take 
negative as well as positive values (one have to keep in 
view that for systems - flows with finite phase velocity 
there is always one zero characteristic Lyapunov expo- 
nent that corresponds to evolution of the infinitesimal 
displacement along the phase trajectory). 

1. If chaos in a Hamiltonian system is suppressed by 
addition of small dissipation, then addition of white noise 
can, as a rule, restore the chaotic motion. The exception 
is found only for very weak chaos, when the system re- 
mains regular at arbitrarily intense noise (see Sec. V.A). 
The increase of the noise intensity raises the Lyapunov 
exponent. In wide context this fact is not quite trivial, 
since there is an example of the system for which noise 
diminishes the (positive) Lyapunov exponent, eventually 
turning it negative |'2l| . 

2. The temperature threshold of chaos depends on 
the strength of chaos of the Hamiltonian system in a 
non-monotonous way (see Fig. 3). For weak chaos it 
increases with the strength (in our model - by linear 
law, see Eq. (|25[) ). since the effective "potential well" 
corresponding to the island of stability becomes deeper, 
whereas for strong chaos the threshold decreases due to 
the shrinking of the islands of stability. 

3. There are two essentially different mechanisms of 
the chaotization of motion by noise. The first one is 
the fluctuational transfer of the motion from the stabil- 
ity island to the locally unstable regions of the phase 
space; its contribution to the Lyapunov exponent de- 
pends on the noise temperature by the " activation law" , 
er+ cx exp(-A/9) (see Eqs. (g3J) and (gSJ and Fig. 4). 
The second is the parametric destabilization inside the 
islands of stability created by small fluctuations of non- 
linear terms of the stroboscopic mapping; its contribution 
to the Lyapunov exponent depends on the noise temper- 
ature by the power law, 0+ cx 6 2 (see Eqs. (|32fl and l|4(J|) 
and Fig. 7). Any one of these mechanisms could be dom- 
inating, depending on the combination of parameters. 

4. Around the threshold of chaos the motion of the 
system with damping and noise differs strongly from the 
chaotic motion of the original Hamiltonian system. It is 
concentrated mainly within the islands of stability with 
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only rare excursions to the domain of the phase space 
occupied by the chaotic component of the prototype. In 
this aspect the restoration of chaos "by noise" differs rad- 
ically from the restoration of chaos "by damping" 
The similarity to the original motion could be reached in 
the domain of strong chaos and high noise temperatures 
6 > 1 (see Fig. 7). 

Lastly it must be noted that the existence of the range 
of parameters where the transition from Hamiltonian to 
dissipative chaos is immediate (see Fig. 1) may be spe- 
cific for the studied model of periodically kicked rotor. In 
this range the influence of noise on chaos has qualitative 



peculiarities: e.g. the increase of noise could reduce the 
Lyapunov exponent (for K — 7 and 7 = 0.1 at 9 = 
a = 1.260(3), and at 9 = 100 a = 1.228(3)). The sce- 
nario of the immediate transition and related problems 
for noisy systems may deserve a special study. 
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